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ABSTRACT 


The  method  of  eliminating  observation  amd  control  spillover 
by  making  groups  of  reduced  order  controlled  modes  orthogonal  to 
each  other  is  investigated.  The  orthogonality  of  these  groups  of 
modes  is  manipulated  by  changing  the  sensor/actuator  model  used  on 
the  structure.  A  sensor/actuator  with  variable  location  and  ori¬ 
entation  is  added  to  the  original  model  to  implement  a  search  for 
a  systematic  method  of  forcing  modes  into  groups. 

Perturbations  in  the  structural  model  are  also  observed  to 
see  how  model  deficiencies  affect  the  stability  of  coupling 
between  modes.  Mass  ch2mges  did  not  affect  the  coupling  signifi¬ 
cantly.  NASTRAN  is  used  to  determine  the  eigenvectors,  which  are 
needed  to  make  the  orthogonality  calculations. 

The  program  ANGLE  is  developed  to  calculate  the  angles 
between  the  modes,  given  a  sensor /actuator  model  emd  the  critical 
modes  taken  from  NASTEIAN.  Groups  of  modes  for  the  decentralized 
controllers  are  identified  and  studied  for  possible  improvements 
that  can  be  achieved  using  sensor/actuator  model  changes. 

After  identifying  likely  Improvable  angles,  the  program 
ORIENT  searches  for  sensor/actuator  models  that  will  indeed  im¬ 
prove  the  desired  angles.  The  ability  to  change  one  angle  without 
significantly  changing  any  others  was  found  to  be  very  difficult. 
The  best  method  for  improving  the  coupling  characteristics  is  an 
iteration  procedure:  improve  one  angle  at  a  time,  then  identify 
the  next  best  way  to  Improve  the  groupings,  and  repeat. 


I  Introduction 


The  ability  of  a  large  space  structure  to  function  as  desired 
is  dependent  on  the  ability  to  control  the  structure.  The  con¬ 
trollability  of  the  structure,  however,  is  hindered  1:^  the  capa¬ 
bility  of  computer  control  systems.  A  solution  to  the  problem  is 
to  use  reduced  order  controllers. 

Reducing  the  number  of  modes  that  are  controlled  introduces 
another  problem.  The  uncontrolled  modes  will  contaminate  sensor 
data.  Diis  contamination  by  the  uncontrolled  modes  produces  ob¬ 
servation  and  ccmtrol  spillover.  Janiszewski  (Ref  1)  shows  that 
the  suppression  of  observation  spillover  allows  the  controller  to 
be  stable.  Spillover  is  eliminated  by  forcing  the  off-diagonal 
terms  of  the  control  matrix  to  equal  zero.  Spillover  terms  can 
also  be  forced  to  zero  by  making  groups  of  modes  orthogonal  to 
each  other.  This  method  of  elimination  was  proposed  by  Miller 
(Ref  2). 

The  modes  of  the  system  matrices  can  be  made  orthogonal  by 
the  selected  placement  of  actuators  and  sensors.  Both  the  orien¬ 
tation  and  location  of  the  actuators  and  sensors  effect  the  ortho- 
gcmality  of  the  modes.  Model  groupings  are  made  of  modes  that  are 
coupled  to  each  other,  but  decoupled  from  the  other  groups  of 
modes.  Miller  suggests  that  the  coupling  characteristics  be  de¬ 
termined  by  calculating  the  angle  between  the  modes  using  a  dot 
product  formula.  Itiis  approach  gives  a  clear  method  of  measuring 
the  coupling  between  modes.  Therefore,  we  have  a  method  for 
determining  the  modal  groupings  ceU-culating  the  coupling  char¬ 
acteristics  due  to  the  sensor/actuator  location  and  orientation 
model. 
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^is  investigation  studied  the  effect  on  the  coupling  of  the 
modes  by  the  cheuiging  of  the  sensor /actuator  model.  The  model  has 
two  parameters  that  are  changed:  the  orientation  and  the  location 
of  the  sensor/actuator  pair.  The  first  method  that  changes  the 
model  is  implemented  by  adding  a  new  sensor/actuator  pair.  This 
method  allows  a  new  and  better  sensor/actuator  model  to  be  chosen 
by  adding  one  sensor/actuator  pair  at  a  time  to  the  model  in  a 
systematic  roeuiner.  This  new  sensor/actuator  can  be  located  at  any 
node  of  the  structure.  Once  a  location  is  chosen,  the  possible 
orientation  angles  can  be  studied.  After  choosing  a  particular 
orientation  (that  gives  favorable  results),  this  sensor/actuator 
becomes  part  of  the  old  model,  and  a  new  sensor/actuator  is  then 
studied.  The  primary  investigation  is  to  choose  a  sensor/actuator 
model  that  will  achieve  an  acceptable  modal  grouping. 

The  second  method  changes  the  sensor/actuator  orientations  of 
the  original  model.  The  determination  of  the  new  orientations 
require  close  observation  of  the  system  matrices.  The  objective 
of  this  investigation  is  to  choose  the  orientations  so  that  the 
existing  modal  grouping  is  improved. 

The  last  area  of  investigation  is  the  effect  the  structural 
modeling  has  on  the  coupling  characteristics  and  modal  groupings. 
The  structure  is  modeled  with  finite  elements.  These  elements 
represent  the  actual  structure.  Thus,  a  structural  change  re¬ 
quires  a  change  in  the  finite  element  model  of  the  structure. 

Changing  the  structural  model  can  vary  the  only  other  para¬ 
meter  used  in  this  study:  the  eigenvectors.  Since  structural 
modeling  is  not  precise,  the  eigenvectors  do  not  accurately  repre¬ 
sent  the  structure.  The  question  is  how  much  effect  does  this 


inaccuracy  have  on  the  control  of  the  structure?  ^e  structural 
analysis  program  NASTRAN  is  used  to  calculate  the  eigenvectors 
using  the  given  finite  element  model.  By  changing  this  finite 
element  model,  inaccuracies  of  the  structural  model  cem  be  tested 
for  their  effect  on  the  coupling  angles  between  modes.  This  is 
done  by  calculating  the  coupling  angles,  using  the  new  eigenvec¬ 
tors  provided  by  the  changed  finite  element  model. 

throughout  this  study,  a  system  model  was  used  with  a  given 
set  of  sensor/actuator  pairs.  Position  sensors  and  force  actua¬ 
tors  are  used.  Thus,  only  translational  degrees  of  freedom  are 
considered  throughout  this  work.  The  sensors  and  actuators  are 
collocated  so  the  modal  coupling  characteristics  will  be  the  same 
for  both  observation  and  control  spillover  terms. 


II  Model 


The  large  space  structure  model  used  in  this  study  is  the 
Charles  Stark  Draper  Laboratory  Model  II  (CSDL-II),  Revision  3. 
The  structure  contains  three  mirrors,  a  focal  plane,  and  the 
suiqport  trusses  for  the  mirrors.  It  also  has  an  equipment  section 
and  two  solar  panels.  The  model  is  shown  in  Fig  1. 

Finite  Element  Model 

Cook  (Ref  3)  defines  finite  element  analysis  as  a  numerical 
procedure  for  solving  a  continuum  mechanics  problem.  Structures 
are  approximated  by  elements  that  are  connected  by  nodes.  The 
CSDL*-II  structure  consists  of  beams  that  form  support  trusses. 
These  beams  are  modeled  by  three-dimensional  frame  elements,  which 
allow  bending,  and  axial  stiffness.  All  of  the  elements  in  this 
model  are  frame  elements. 

The  elements  in  the  model  represent  hollow  graphite  epoxy 
tubes.  The  wall  thicknesses  of  the  tubes  range  from  0.03cm 
0.067cm.  The  radius  ranges  from  3.6cm  to  8.1cm.  Each  element  is 
designed  so  that  it  is  only  strong  enough  to  meet  local  buckling 
and  member  natural  frequency  constraints.  The  natural  frequency 
of  each  element  has  to  be  greater  than  lOhz  so  that  local  vibra¬ 
tions  will  not  interact  with  system  vibrations.  For  those  inter¬ 
ested,  the  buckling  and  natural  frequency  constraints,  and  the 
sizes  and  section  properties  of  each  element  are  given  by  Hender¬ 
son  (Ref  4:3-6);  however,  these  properties  were  not  used  in  this 
investigation. 


4 


5 


Each  element  is  connected  to  another  by  a  single  node.  The 
structure  has  a  total  of  59  nodes  with  23  lumped-mass  nodes. 
These  nodes,  along  with  their  masses,  are  cataloged  in  Table  I. 
The  locations  of  the  nodes  are  shown  in  Fig  2. 

The  mass  model  is  a  revision  of  the  previous  model  used  for 
CSDL-II.  In  this  revision,  the  mirror  masses  of  the  model  are 
changed  so  that  details  of  the  support  truss  of  the  rigid  mirrors 
can  be  modeled  (Ref  4).  The  total  mirror  mass  of  the  previous 
model  is  kept  the  same.  However,  the  toted  mass  of  each  mirror  is 
now  divided  into  two  parts,  the  mirror  mass  and  the  optical  sup¬ 
port  structure  mass.  Thus,  nodes  1001,  1002,  and  1003,  the  nodes 
at  the  center  of  each  mirror,  contain  the  lumped  mirror  masses. 
The  nodes  surrounding  the  mirrors  have  the  lumped  masses  due  to 
the  optical  support  structure  mass.  Both  the  primary  and  tertiary 
mirrors  already  have  a  massive  support  truss  at  one  end.  Thus, 
this  support  mass  has  already  been  partially  accounted  for. 
Therefore,  the  mass  of  the  support  element  is  subtracted  from  the 
lumped  optical  support  mass.  This  leads  to  large  lumped  masses 
(due  to  the  optical  support  mass)  at  two  corners  and  smaller 
masses  (due  to  the  optical  mass  minus  the  support  truss  mass)  at 
the  other  two  corners  of  each  mirror,  since  that  is  where  the 
massive  support  truss  is  located. 

The  finite  element  analysis  is  accomplished  through  the 
structural  analysis  program  NASTRAN.  Each  element  has  stiffness 
properties,  and  some  nodes  have  lumped-mass  properties.  Mass  and 
stiffness  matrices  for  the  system  are  calculated  by  NASTRAN. 
These  matrices  allow  NASTRAN  to  complete  an  eigenvector  analysis. 
These  eigenvectors  are  then  used  to  calculate  the  coupling  charac¬ 
teristics  of  the  structure. 


nie  NASTSAN  bulk  data  deck  Is  listed  in  J^pendix  A.  Ttie  bulk 
data  deck  produces  a  listing  of  the  eigenvalues  and  eigenvectors 
for  the  finite  element  model.  To  change  the  modeling  of  the 
structure,  the  data  deck  is  modified.  New  eigenvalues  2md  eigen¬ 
vectors  will  be  produced  with  each  change  in  the  data  deck. 

The  mass  model  of  CSDL-II  is  modified  by  changing  the  lumped 
mass  values  in  the  data  deck.  This  is  accomplished  through  the 
CONM-2  statements.  There  is  one  CONH-2  line  for  each  node  that 
has  a  lumped  mass  (Appendix  A).  Bach  line  contains  the  node 
number  of  the  lumped  mass  and  the  V2due  of  the  mass. 

Sensor /Actuator  HfijSel 

The  sensor/actuator  pairs'  location  and  orientation  are  given 
in  Table  II.  This  model  was  suggested  by  CSDL  (Ref  4).  The  sug¬ 
gestion  is  based  on  the  line-of-sight  (LOS)  error.  The  LOS  error 
is  derived  from  the  mathematical  model  given  by  CSDL  (Ref  5).  Ihe 
error  is  a  function  of  21  coordinate  variables.  These  21  vari¬ 
ables  are  the  coordinates  of  12  nodes  that  control  the  geometry  of 
the  optical  path.  Thus,  to  control  the  LOS  error,  these  21  coor¬ 
dinates  were  chosen  for  the  sensor/actuator  pairs  locations. 
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laor /Actuator  Nodes  and  Orient 


ALPHA  BETA 

90.00  0.00 

90.00  90.00 

90.00  90.00  0.00 

0.00  90.00  90.00 

90.00  0.00  90.00 

90.00  90.00  0.00 

90.00  90.00  0.00 

0.00  90.00  90.00 

90.00  0.00  90.00 

90.00  90.00  0.00 

90.00  90.00  0.00 

90.00  0.00  90.00 

90.00  90.00  0.00 

90.00  90.00  0.00 

90.00  90.00  0.00 

90.00  90.00  0.00 

0.00  90.00  90.00 

90.00  0.00  90.00 

90.00  90.00  0.00 

90.00  0.00  90.00 

90.00  90.00  0.00 
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Note  that  all  of  the  sensor /actuator  pairs  are  oriented  along 
2m  axis  direction  as  specified  by  the  LOS  equation.  This  restric¬ 
tion  is  not  followed  with  sensor/actuator  pairs  that  are  added 


I 

I 

throughout  the  study. 
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Ill  Theory 


Equations  Motion 

A  vibrating  structure  has  equations  of  motion 

M^+Eg  +  Kg  =  Q  (1) 

where  g  is  a  n-vector  of  generalized  coordinates,  H  a  n  x  n  sym¬ 
metric  mass  matrix,  K  a  n  x  n  symmetric  stiffness  matrix,  E  a  n  x 
n  symmetric  damping  matrix,  and  Q  a  n-vector  of  input  forces. 
'Diis  equation  is  also  valid  for  large  space  structures  such  as  the 
CSDL-II  model. 

The  input  force  vector  can  be  separated  into  two  parts 

Q  =  D  u  (2) 

where  D  is  a  n  x  m  matrix  of  direction  cosines,  and  u  is  a  m- 
vector  of  actuator  force  inputs.  The  direction  cosine  matrix,  D, 
is  block  diagonal  when  there  is  only  one  actuator  per  node.  The 
matrix  columns  contain  the  three  direction  cosines  of  the  corres¬ 
ponding  actuator  force  inputs.  Equation  1  is  simplified  by  using 
modal  coordinates,/^.  Janiszewski  shows  that  the  equations  can  be 
expressed  in  state  variable  form 


where 


the  general  output  equation  is 


y  =  C  X 


(4) 


where 

C  “  I  ] 

The  modal  matrix/  ^  /  Is  a  matrix  containing  the  eigenvectors 
of  the  structure  by  column.  It  is  a  matrix  which  is  n  x  n^  where  n 
is  the  number  of  degrees  of  freedom  and  n^  the  number  of  controlled 
eigenvectors. 

Since  we  only  have  point  sensors  and  force  actuators  that  are 
collocated/  half  of  each  matrix  is  filled  with  zeros.  If  the  zero 
halves  of  these  matrices  are  neglected/  we  may  define 


and  ng  x  n^respectively,  where  is  the  number  of  actuators  and 
ng  is  the  number  of  sensors. 


r# 


Coupling 

In  a  complex  structure  where  there  are  a  large  number  of 
modes,  a  reduced  order  model  is  used.  Janiszewski  shows  that  only 
a  subset  of  the  total  number  of  structural  modes  need  to  be  con¬ 
trolled.  If  this  subset  of  modes  is  still  large,  then  multiple 
controllers  are  used. 

Each  controller  is  designed  to  operate  on  only  a  small  number 
of  the  modeled  modes,  which  are  a  subset  of  the  structural  modes. 
However,  spillover  terms  occur  because  the  controller  actually 
operates  on  all  the  modes.  l!hus,  these  spillover  terms  need  to  be 
eliminated. 


Bi  Gj  =  0 


i^  j 


(6) 


where  G  is  the  control  feedback  gain  matrix,  and  K  is  the  observer 
gain  matrix.  Hie  subscripts  on  the  matrices  correspond  to  groups 
of  modes. 

Hiller  shows  that  equations  6  are  satisfied  if  the  system 
matrices,  B  and  C  of  each  controller,  can  be  made  orthogonal  to 
the  remaining  controllers.  Therefore,  the  spillover  terms  are 
eliminated  if  groups  of  modes  are  found  with  the  following  pro¬ 
perty:  All  modes  in  a  group  are  coupled,  and  each  group  is  decou¬ 
pled  from  the  other  groups. 


16 


The  modes  that  have  been  referred  to  are  actually  the  rows  of 
B  and  the  columns  of  C.  Therefore,  it  is  the  rows  of  B  (columns 
of  C)  that  must  be  separated  into  modal  groups.  These  rows 
(columns)  are  the  eigenvectors  multiplied  by  the  direction  cosine 
matrix. 

The  degree  of  coupling  is  measured  by  the  angle  between  the 
modal  vectors.  For  B,  the  angles  between  the  rows  are 


S  »  COS* 


BL  ■  Bj 

IIbJI  Iliill 


(7) 


where  and  Bj  are  rows  i  and  j  respectively  of  B  and  ||  B^jlis  the 
norm  of  the  1^^  row  vector  of  B.  These  angles  can  be  made  into  a 
symmetric  n^  x  n^  matrix  of  angles. 

For  C,  the  angle  is  taken  between  the  columns.  With  the 
collocation  of  sensors  and  actuators,  the  angle  matrix  produced 
from  B  amd  C  are  identical. 


Location 

Ihe  elements  of  B  and  C  are  functions  of  the  location  of  the 
actuators  and  sensors.  An  element  is  determined  from  three  values 
of  an  eigenvector  from  the  modal  matrix  ^  which  are  multiplied  by 
the  direction  cosines  of  an  actuator  or  sensor.  Therefore,  an 
element  of  B  is 

®ij  “  (^6  •  node(j)-5  i  COSCXj  +  .  node(j)-4  i  (8a) 

+  •  node(j)-3  i  COSjTj 

where  node  (j)  is  the  node  of  the  j'th  actuator.  If  k  »  6» 
node(j)-5,  then  the  number, is  the  k'th  row  element  of  the 
i'th  eigenvector ,  in  the  modal  matrix,  0.  These  matrices  can 
also  be  determined  as  follows: 
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where  k  »  1,  13/  7  for  the  first  three  actuators  in  this  example 
In  this  case/  the  three  actuators  are  at  nodes  1/  3,  and  2,  respec 
tively.  ThuS/  k  for  actuator  2  is  calculated  by  6x(3)-5. 


Similarly r  an  element  of  C  is  given  by 

^ij  “06  •  node(i)-5  j  <»SO^i  +  (j>^  .  node(i)-4  j  ^^^i 
06  •  node (i) -3  j  ^^i 

Note  that  by  adding  an  actuator ^  a  column  is  added  to  B. 
With  a  new  sensor,  a  row  is  added  to  C.  Also,  if  an  extra  eigen¬ 
vector  is  added  to  the  modal  matrix  (ji ,  then  a  row  is  added  to  B 
and  a  column  to  C. 

The  magnitude  of  the  elements  of  B  and  C  is  dependent  on  the 
magnitude  of  the  eigenvectors  with  respect  to  the  norm.  If  the 
eigenvectors  have  small  values  with  respect  to  the  norm  at  the 
location  of  the  j'th  actuator  and  i'th  sensor,  then  the  elements 
of  B  and  C  corresponding  to  the  actuator  and  sensor  will  be  small. 
So  an  element  of  the  system  matrices,  B  and  C,  can  be  made  small 
if  the  actuator/sensor  location  is  chosen  so  that  the  eigenvector 
value  is  small  with  respect  to  its  norm.  Likewise,  if  the  value 
of  the  eigenvector  is  large  with  respect  to  its  norm  at  a  parti¬ 
cular  location,  then  the  system  matrix  element  will  be  large. 

The  angle  between  modes  is  dependent  on  the  elements  of  the 
system  matrices.  Therefore,  the  angle  is  also  a  function  of  the 
location  of  the  sensor/actuator  pairs.  Desired  angles  can  be 
achieved  by  the  judicious  placement  of  sensor/actuator  pairs.  The 
addition  of  a  sensor/actuator  pair  will  allow  the  angles  between 
modes  to  be  varied,  since  a  new  emd  different  element  is  added  to 
each  mode  of  the  system  matrices.  Thus,  the  angle  between  modes 
is  dependent  on  the  system  matrix  which  includes  the  new  elements 
which  are  dependent  on  the  new  sensor  actuator  pair. 
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Adding  an  element  to  each  mode  of  the  system  matrices  will 
yield  the  following  equation  for  the  angles  between  two  modes 
(using  equation  7  with  A  and  B  the  two  modes): 


0  =  CX>S 


rf 


/\i  Bi  *  AtBi  ♦ . . .  ♦  Ak  Bx 


(10) 


%diere  A^  and  are  the  new  elements  of  each  mode. 
7t:is  equation  can  be  simplified 


S^cos‘ 


•  B  Ax  Bx 


(11) 


where 


A‘  =  Z4‘  &"  =  'Zb\ 

i*l  i*l 

the  dot  product,  A  •  B  is  taken  without  the  new  elements. 

Equation  (11)  can  be  reduced  if  one  of  the  new  elements  is 
equal  to  zero.  Letting  B^  equal  zero,  we  have 


0=  COS  -  II  g  II  J 


(12) 


If  A^  is  also  equal  to  zero,  the  angle  does  not  change. 
However,  the  angle  will  change  if  A,^  is  not  zero.  No  matter  what 
value  A^  is  (the  larger  the  better),  the  denominator  of  equation 
(12)  must  increase.  Thus,  this  leads  to  the  angle  going  to  90°. 
Therefore,  to  decouple  two  coupled  modes,  find  a  location  where 
one  mode  has  large  (relative  to  norm)  values  and  the  other  small 
values. 
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Olie  cheuige  in  the  cosine  of  the  angle  can  be  determined  by 
setting  the  previous  cosine  equal  to  the  new  cosine  of  the  angle 
(equation  12) . 


COS0= 


”  A-  B 

A  '  B 

Pll  IIBIlJ 

Ulus,  the  cosine  of  the  angle  will  change  the  factor 

II  All  _  Wf 

A 

“  (A'tAi)'-^ 


If  is  small  compared  to  the  norm  (taken  without  A^),  the 
cosine  of  the  angle  will  not  change. 

A  second  case  occurs  if  the  new  element  product  dominates  the 
previous  dot  product.  This  case  can  occur  when  the  original  modes 
are  decoupled r  since  the  dot  product  is  zero.  The  angle  now  re¬ 
duces  to  —  _ 


0  =  cos 


A,  8. 


(A"*A.*r(B'*  Bi)'' 


r 

1 

1  1 

r' 

Vi 

[Ki)'] 

fit) 


The  new  angle  approaches  zero  if  the  new  elements  are  large 
compared  to  the  previous  norms.  Thus,  two  modes  can  be  coupled  if 
a  location  for  the  sensor/actuator  pair  can  be  found,  such  that 
the  new  elements  are  large  with  regard  to  the  norm.  Therefore, 


one  way  of  coupling  two  modes  is  to  put  a  sensor /actuator  pair  at 


a  node  where  the  eigenvector  is  large.  Also,  if  the  elements 
produced  for  the  other  modes  are  small,  then  only  the  modes  for 
which  the  elements  are  large  will  have  a  significant  change  in 
their  coupling  characteristics. 

There  are  two  disadvantages  to  this  method.  First,  even 
though  the  emgles  between  two  modes  with  two  new  small  elements 
remain  the  same,  the  angle  between  modes  with  new  small  elements 
and  modes  with  new  large  elements  drastically  change  (they  decou¬ 
ple  as  shown  in  equation  (12)).  Secondly,  there  is  no  guarantee 
that  a  node  location  and  orientation  will  be  found  that  will 
satisfy  these  angle  c<mdltlons.  Hiis  is  because  the  location  and 
orientation  simply  determine  which  values  within  the  eigenvector 
are  used.  If  the  eigenvectors  do  not  contain  precisely  the  right 
values,  it  is  impossible  to  pick  a  location  and  orientation. 

€• 

Oiientation 

The  location  of  the  sensor/actuator  pair  determines  which 
elements  (determined  by  k)  of  the  eigenvector  are  used  in  the 
element  equation  (8).  Die  three  elements  chosen  by  the  node  loca¬ 
tion  of  the  sensor/actuator  are  taken  from  each  eigenvector  and 
multiplied  to  the  direction  cosines  of  the  sensor /actuator.  After 
this  calculation  the  elements  of  B  then  become  only  a  function  of 
the  orientation  (through  matrix  D).  Thus,  the  set  of  orientation 
angles  must  be  known  so  the  element  can  be  determined. 

The  orientation  angles,  ,  ^3  r  and  }S  are  related.  If  any 
two  and  the  octant  of  the  sensor/actuator  are  known,  then  the 
third  can  be  determined.  If  o<  and  0  are  known,  ^  is  then  deter¬ 
mined  as  shown  in  Fig  3. 
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cos  cx  =  OA/r  = 
cos  /3  =  OB/^  =  y  (15) 

cos  =  6c/r  =  z 

since  P  is  a  unit  vector,  the  direction  cosines  are  equal  to 
the  components  of  the  vector  along  the  axes. 

Die  magnitude  of  P  is: 

Hr  II  =  -Ix’.y**  z‘  ‘  =  I  (16) 

Substituting  equations  (15)  into  equation  (16) 

COs2o(+  00s26+  cos  1  (17) 

Solving  for  5  #  we  have 

^  =  COS”l  J  1  -  Icos2o<+  cos2^]  '  (18) 

where  0(  <md  are  varied  from  0^  to  180^. 

All  possible  orientation  angles  can  be  tried,  or  specific 
orientations  can  be  selected,  so  that  existing  elements  of  the  B 
and  C  matrices  are  changed.  The  challenge  in  this  procedure  is  to 
decide  what  the  elements  of  these  matrices  should  become  in  order 
for  the  modal  angles  to  be  acceptable.  Dien  specific  orientations 
are  chosen  so  that  the  desired  elements  of  the  system  matrices  are 
produced. 

Decoupled  modes  have  a  dot  product  equal  to  zero  (relative  to 
their  norms).  Thus,  every  large  element  of  one  mode  will  have  a 
corresponding  term  in  the  other  mode  that  is  small  (relative  to 
norms).  This  correspondence  between  the  modes  causes  the  dot 
product  to  go  to  zero,  which  leads  to  decoupled  modes.  Another 
possibility  is  that  large  products  have  negative  products  that 
cancel  them. 
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Coupled  inodes  have  the  opposite  condition:  large  elements 
correspond  to  each  other r  and  small  elements  correspond  to  each 
other.  fHiuSr  the  dot  product  equals  the  product  of  the  norms. 

By  changing  the  sensor /actuator  orientationr  these  conditions 
can  be  switched.  The  mode  elements  for  the  decoupled  modes  need 
to  cause  the  dot  product  to  increase  in  magnitude  or  the  norms  to 
decrease.  This  is  accomplished  by  causing  an  element  pair  (one 
corresponding  large  and  small  element)  to  become  both  small  or 
both  large. 

The  coupled  modes  are  decoupled  by  decreasing  the  dot  product 
or  increasing  the  norms.  Again,  this  is  accomplished  by  causing 
an  element  pair,  that  has  either  both  large  or  both  small  elements 
to  readjust,  so  that  the  element  pair  has  one  large  element  and 
one  small  element.  Note  again,  all  magnitudes  above  are  relative 
to  the  norms. 

Modes 

Twelve  eigenvectors  were  suggested  1:^  LHSC  (Ref  6)  to  be  the 
critical  modes.  These  modes  are  listed  in  Table  III.  The  modes 
were  selected  according  to  three  criteria:  modal  cost,  control¬ 
lability,  and  observability.  Ten  modes  were  selected  that  had  the 
highest  controllability/observability  product.  Also,  two  modes 
were  selected  with  the  highest  modal  cost.  These  modes  were  used 
throughout  this  stud|y  as  the  critical  modes. 

These  modes  can  change  when  the  eigenvectors  themselves 
change.  With  each  perturbation  of  the  finite  element  model,  the 
eigenvectors  (^matrix)  may  change.  'Aiis  change  obviously  changes 
the  modes  (B  and  C  matrices).  This  was  the  only  change  to  the 
modes  that  was  studied. 


Table-Ill 


node 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 
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coaputer  programfi 

The  computer  programs  used  in  this  study  are:  NASTRAN, 
SELECT,  ANGLE,  and  two  versions  of  ORIENT.  As  stated  previously, 
NASTRAN  is  used  to  perform  a  finite  element  analysis  of  CSDL-II. 
The  program  outputs  the  eigenvectors  and  eigenvalues  into  a  per¬ 
manent  file,  flhis  file  of  eigenvectors  is  the  basis  for  the  other 
programs. 

SELECT  is  listed  in  ^pendix  b.  The  purpose  of  this  program 
is  to  pick  off  the  eigenvectors  output  by  NASTRAN.  Any 
combination  of  eigenvectors  in  any  order  can  be  selected.  The 
eigenvectors'  numbers  are  first  read  by  the  program.  The  file  of 
eigenvectors  is  stored  in  a  three-dimensional  matrix.  Then,  the 
selected  eigenvectors  are  put  into  an  output  permanent  file. 

ANGLE  is  listed  in  Appendix  C.  ANGLE  inputs  the  modal  vector 
(from  SELECT),  and  a  sensor/actuator  model  that  has  been  chosen. 
It  then  calculates  the  B,  C,  and  angle  matrices.  The  B  and  C 
matrices  are  calculated  from  the  element  equations  (8-9).  The 
angle  matrix  is  determined  by  using  the  cosine  equation  (7).  The 
angle  matrix  is  then  used  to  determine  modal  groupings. 


ORIENT  is  the  program  that  determines  the  orientation  for 
each  node  location  that  will  meet  the  desired  requirements.  Two 
diverse  versions  of  ORIENT  are  used  to  determine  the  acceptable 
orientations.  Both  work  by  eliminating  those  orientations  that  do 
not  meet  the  criteria  (below).  Version  two  is  listed  in  i^pendix 
D. 

The  first  version  of  ORIENT  tries  to  find  the  orientations 
that  will  cause  large  (relative  to  the  norms)  elemental  values  in 
two  modes,  but  small  elemental  values  in  all  other  modes.  It  does 
this  by  checking  the  values  of  each  new  element.  If  the  two  modes 
in  question  were  previously  coupled,  they  would  then  be  decoupled. 

The  second  version  simply  prints  all  orientations  that  have 
the  modes  in  a  given  set  of  modal  groupings.  The  criteria  used, 
if  the  modes  are  to  be  coupled,  is  that  the  angle  must  be  less 
than  82®  or  greater  than  103®,  If  the  modes  are  to  be  decoupled, 
then  the  angle  must  be  between  86®  and  97®.  Obviously  it  is  de¬ 
sired  to  have  coupled  angles  at  0®  or  180®,  and  decoupled  angles 
at  90®,  but  it  is  not  possible  to  move  all  the  modal  angles  to 
these  desired  angles.  These  specific  angles  were  chosen  so  that 
already  acceptable  modal  angles  would  fall  within  their  limits  and 
not  have  to  change.  Yet  the  worst  angles  would  have  to  change. 
Note  that  this  program  is  very  versatile  because  conditions  can  be 
(hanged  so  that  almost  any  objective  can  be  investigated.  If  the 
modal  groupings  need  to  be  changed,  all  that  has  to  be  changed  is 
the  input  mode  numbers.  Similarly,  the  bounds  for  coupling  and 
decoupling  (»n  be  easily  chemged. 


IV  Investigation 


SaiBOg/Actuatgr  tlodel 

The  thrust  of  the  investigation  is  to  observe  the  coupling 
characteristics  of  the  B  and  C  matrices  due  to  changes  in  the 
sensor /actuator  model.  Thus,  determining  matrices  B  and  C  is  the 
first  objective.  Using  the  given  sensor/actuator  model,  the 
matrices  are  calculated. 

These  matrices  are  found  through  the  ANGLE  program.  The 
program  was  run  originally  with  17  modes;  however,  five  eigenvec¬ 
tors  were  dropped:  11,  14,  16,  26,  29  (modes  13,  14,  15,  16,  17). 
This  was  done  to  simplify  the  problem  (fewer  modal  angles).  The 
five  modes  dropped  had  the  lowest  modal  cost  and  controllability/ 
observability  product  of  the  17  eigenvectors.  Thus,  the  12  modes 
left  were  used  as  critical  modes  throughout  the  rest  of  the  inves¬ 
tigation  . 

Next,  the  angle  between  the  modes  was  calculated  using  ANGLE. 
The  results  are  in  the  form  of  a  symmetric  matrix  where  the  rows 
and  columns  correspond  to  the  mode  numbers.  Therefore,  a  12  x  12 
angle  matrix  is  determined  (see  Table  IV). 

The  purpose  of  calculating  the  angles  between  the  modes  is  to 
determine  modal  groupings.  These  groupings  are  then  used  by  the 
different  controllers.  The  original  sensor/actuator  model  pro¬ 
duced  three  groups  of  modes  that  are  in  general  orthogonal  to  each 
other  or  coupled  within;  however,  the  angles  do  not  necessarily  do 
both.  Since  modes  11  and  12  are  decoupled,  a  second  modal  group¬ 
ing  is  also  possible.  These  two  possible  groupings  are  shown  in 


Table  V 


Both  groupings  have  problems.  The  problem  angles  of  the 
first  group  are  circled  in  Table  IV,  while  the  second  group's 
problem  angles  are  boxed.  Hheae  angles  are  a  problem  because  they 
are  not  as  nicely  coupled  or  decoupled  as  the  other  angles.  The 
first  grouping  has  decoupled  modes  11  and  12  in  group  II,  where 
they  should  be  coupled.  The  angles  between  modes  2-5,  2-9,  4-5, 
4-9,  and  8-10  indicate  coupling,  but  these  modes  are  in  different 
groups  (so,  should  be  decoupled).  The  second  grouping  puts  mode 
12  into  group  III.  While  this  separates  modes  11  and  12,  mode  12 
is  still  coupled  to  modes  1,  6  and  7  in  group  II. 
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Table  IV 


Angle  Matrix  with  Original  Sensor /Actuator  Model 


1 

0.0000 

89.9876 

89.9991 

2 

89.9876 

0.0000 

67.4653 

3 

89.9991 

67.4653 

0.0000 

4 

90.0087 

175.3876 

115.3141 

5 
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fll6.852j 

V  i 

1  104.0098 

6  123.2047 

90.0533 

89.9855 

7  107.4792 

89.9662 

90.0376 

8 

90.0478 
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116.9917 

9 

90.0128 

^16.232| 

1  91.5675 

.0 

90.0066 

90.1561 

80.1261 

.1 

39.1619 

87.8879 

89.4097 

.2 

81.1169 

87.9411 

89.2693 
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S. 

1 

A 

90.9987 


5.  A 

90.0884  123.2047 
90.0533 


115.3141 

104.0098 

89.9855 

0.000 

[60.6688 

1  89.9526 

[60.6688 

1  0.0000 

89.6936 

89.9526 

89.6936 

0.0000 

90.0220 

90.0156 

61.2182 

73.1913 

102.7661 

90.0720 

66.4898' 

1  68.5695 

89.8722 

90.6814 

111.8732 

90.1536 

92.4490 

96.6087 

104.9007 

92.3517 

95.6278 

73.2904 

10. 

11 

12 

1  107.4792 

2  89.9662 

3  90.0376 

4  90.0220 

5  90.0156 

6  61.2182 

7  0.0000 

8  89.9103 

9  90.1494 

.0  89.9571 

.1 
.2 


80.5613 

tl27.804l| 


90.0478 

90.0128 
- - 

90.0066 

106.5687 

^16.232^ 

90.1561 

116.9917 

91.5675 

80.1261 

73.1913 

(66.4898] 

1  90.6814 

102.7661 

68.5695 

111.8732 

90.0720 

89.8722 

90.1536 

89.9103 

90.1494 

89.9571 

0.0000 

93.8830 

(38.0463] 

93.8830 

0.0000 

103.2217 

[38.0463} 

103.2217 

0.0000 

93.4985 

90.5118 

94.1486 

93.1022 

90.7436 

93.4951 

39.1649 

81.1169 

87.8879 

87.9411 

89.4097 

89.2693 

92.4490 

92.3517 

96.6087 

95.6278 

104.9007 

73.2804 

80.5616  [ 

127.804]| 

93.4985 

93.1822 

90.5112 

90.7436 

94.1486 

93.4051 

0.0000 

(mIsiw 

89.8169 


0.0000 


The  two  groupings  represent  two  approaches  to  the  problem  of 
modes  11  and  12.  The  first  grouping  is  an  attempt  to  couple  the 
two  modes  together.  The  second  grouping  places  mode  12  into  group 
III.  This  grouping  tries  to  decouple  mode  12  from  the  other  modes 
in  group  II.  This  first  approach  is  the  simpler  of  the  two.  It 
is  easier  to  couple  modes  together,  since  all  that  is  needed  is  a 
change  of  about  10°  off  90°.  Although  this  is  not  much  of  a 
change,  it  is  probably  good  enough  for  the  controller  to  work. 

To  couple  the  modes  together,  the  sensor/actuator  model  is 
changed.  The  original  model  stayed  constant,  but  a  new  sensor/ 
actuator  pair  was  added.  This  pair  could  be  placed  at  any  loca¬ 
tion  and  at  any  orientation.  The  decision  on  where  to  place  the 
sensor /actuator  and  at  what  orientation  is  complicated.  This  com¬ 
plication  arises  because  every  location  and  orientation  chosen 
will  change  every  angle  between  the  modes. 

To  combat  this  problem,  the  two  versions  of  ORIENT  pick  out 
the  orientations  that  will  couple  modes  11  and  12,  but  leave  the 
others  in  the  same  modal  groupings. 

Both  versions  try  every  node  point  location.  This  is  because 
there  is  no  way  of  knowing  which  node  point  is  the  best.  Attempts 
to  rank  the  nodes  were  not  precise  enough  to  eliminate  the  other 
nodes  from  consideration.  Therefore,  it  was  better  to  look  at 
every  node  than  to  discard  any  of  the  nodes. 

After  the  program  picks  a  node,  ORIENT  checks  all  possible 
orientations  which  will  satisfy  the  coupling  conditions.  ORIENT 
checks  the  orientations  by  calculating  the  modal  angles  until  one 
is  rejected.  If  an  angle  is  rejected,  the  next  orientation  is 


tried;  if  all  angles  are  accepted,  the  emgle  matrix  is  then  prin¬ 
ted.  Host  of  the  orientations  are  discarded,  but  the  output  still 
has  many  orientations  that  will  satisfy  minimum  conditions.  T!tie 
output  orientations  from  ORIENT  were  then  checked  by  hand  to 
determine  the  best  orientation  that  gives  acceptable  groupings. 
The  determination  of  the  best  orientation  necessarily  chooses  the 
best  node  location. 

The  next  consideration  of  sensor/actuator  model  changes  is 
that  of  changing  the  grouping.  The  first  grouping  was  chosen 
because  it  was  an  acceptable  grouping  of  the  original  sensor/ 
actuator  model.  This  grouping  and  the  second  one  described  above 
are  not  the  only  acceptable  groupings.  It  does  not  matter  which 
modes  are  in  a  given  group,  as  long  as  they  are  coupled  to  each 
other,  and  decoupled  from  all  other  modes.  This  leads  to  an 
investigation  of  what  other  modal  groupings  are  possible  with  a 
sensor/actuator  model  change. 

Again,  modes  11  and  12  were  considered.  ORIENT  was  run  with 
two  cases.  Case  one  is  that  modes  11  and  12  be  decoupled.  Case 
two  is  that  modes  11  and  12  be  coupled.  No  other  restrictions  are 
made.  This  gives  an  output  of  all  possible  modal  grouping  where 
11  and  12  are  either  coupled  or  decoupled. 

Ihe  last  type  of  model  change  investigated  was  changing  the 
original  sensor /actuator  orientations.  This  modal  change  was  in¬ 
vestigated  to  see  if  the  original  orientations  are  indeed  the 
best.  Two  extremes  were  studied.  First,  change  the  orientation 
so  that  decoupled  modes  are  coupled.  Second,  decouple  the  coupled 
modes.  Ihe  hope  was  that  only  the  angles  of  the  modes  in  question 
would  be  radically  changed.  Ohe  modal  emgles  (old),  the  coupling 
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objectives,  cmd  the  new  sensor/  actuator  orientations  are  listed 
in  Table  VI.  In  each  run  of  the  ANGLE  program,  only  one  sensor/ 
actuator  orientation  was  changed.  Again,  this  was  done  to  keep 
other  angles  from  changing  too  much. 
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TABLE  VI 


I  Original  Sensoc/Afituator  Model  Orientation  Changes 


MODES 

OLD 

ANGLE 

OBJECT 

ACTCMOB 

ORIENTATION 
ALPHA  BETA 

CHANGES 

GAMMA 

11-12 

89.8° 

Couple 

2 

0° 

90° 

90° 

11-12 

89.8° 

Couple 

13 

90° 

45° 

45° 

11-12 

89.8° 

Couple 

10 

90° 

45° 

45° 

8-10 

38° 

Decouple 

5 

45° 

45° 

90° 

8-10 

38° 

Decouple 

10 

45° 

90° 

45° 

8-10 

o 

00 

Decouple 

12 

90° 

45° 

45° 

8-9 

93.8° 

Couple 

21 

45° 

90° 

45° 

7-10 

89.4° 

Couple 

4 

45° 

45° 

90° 

7-10 

89.4° 

Couple 

3 

90° 

0° 

90° 
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MasB  Model  Chan9eB 


Determining  the  sensor/actuator  pairs' location  and  orien- 
I  tation  that  produce  a  good  grouping  of  modes  is  only  valid  for  the 

structural  model  that  has  those  eigenvectors  used  in  the  deter¬ 
mination.  If  the  structural  model  is  in  error,  then  the  eigen- 
I  vectors  are  in  error.  This  error  in  the  eigenvectors  may  change 

the  coupling  angles  between  the  modes.  If  the  coupling  character¬ 
istics  are  drastically  changed,  the  modal  grouping  could  be  dif- 
i  ferent.  This  change  in  the  grouping  is  unacceptable. 

This  problem  led  to  finite  element  model  changes  being  stu¬ 
died.  Five  different  model  ch^ulges  were  made.  These  are  shown  in 

* 

Table  VII.  The  first  three  are  error  percentage  changes  in  the 
total  mass  of  the  tertiary  mirror.  A  percentage  of  the  total  mass 
is  added  in  proportion  to  each  node.  The  fourth  model  keeps  the 
I  a  same  total  mass  but  changes  the  proportion  of  mass  at  each  node. 

The  fifth  model  adds  mass  to  the  nodes  with  smaller  mass  amounts. 
The  mass  amount  added  equalizes  the  mass  at  each  corner  node, 
i  All  five  mass  changes  were  on  the  tertiary  mirror.  This 

mirror  was  chosen  so  that  the  structural  model  could  be  tested 
without  symmetry  effects.  Since  only  this  mirror  would  receive 
mass  changes,  no  symmetry  sensitive  eigenvectors  would  go  unaffec¬ 
ted  by  the  changes.  Also,  the  three  different  types  of  mass 
changes  further  decreased  this  possibility. 

NASTRAN  was  then  run  with  each  of  the  new  models.  Ihis  pro¬ 
duced  new  eigenvectors  for  each  model.  After  picking  off  the 
critical  eigenvectors,  ANGLE  was  run  to  determine  the  coupling 
angles  between  the  modes.  The  sensitivity  of  the  modal  groupings 
to  the  mass  modeling  changes  made  to  the  tertiary  mirror  was  then 
observed. 
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V  ResUltB 


The  results  of  this  investigation  come  from  two  sources. 
First,  the  sensor/actuator  model  is  manipulated  with  the  objective 
of  obtain  an  acceptable  grouping  of  modes.  Further,  the  model  is 
chemged  again  to  improve  the  selected  groupings.  Also,  changes  of 
the  original  orientations  are  studied.  Second,  the  finite  element 
model  is  perturbed  so  mass  modeling  errors  can  be  observed. 
Specifically,  the  lumped  masses  at  the  nodes  are  changed. 

Sensor /Actuator  tfodel 

The  original  sensor/actuator  model  is  listed  in  Table  II. 
The  ANGLE  program  was  run  using  this  data.  The  angles  between  the 
modes  are  given  in  Table  IV.  By  observing  these  coupling  angles, 
modal  groupings  were  picked  out.  The  two  groupings  chosen  are 
listed  in  Table  V. 

Ilie  results  of  ORIENT  runs  are  mixed,  when  requirements  of 
acceptance  (see  Appendix  D)  were  too  strict,  no  sensor/actuator 
models  were  found.  Ihls  predicament  was  caused  by  some  modal  cou¬ 
pling  angles  being  too  far  off  the  desired  angles.  In  other 
words,  the  angles  were  required  to  be  so  close  to  either  coupled 
(0^  or  180^)  or  decoupled  (90^)  that  some  bad  angles  could  not 
move  inside  the  limits.  These  are  the  problem  angles  listed  in 
the  investigation  section.  Thus,  these  particular  problem  angles 
were  then  excluded  from  the  acceptance  requirements  of  ORIENT. 
The  only  exception  is  the  angle  between  modes  11  and  12.  Since 
this  angle  was  one  of  the  worst,  it  was  singled  out  to  continue  to 
try  to  meet  the  acceptance  limits.  Ihe  results  from  ORIENT  with 
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these  new  requirements  of  acceptance  (most  problem  angles  ex¬ 
cluded)  give  node  42  as  the  best  location.  There  were  fifteen 
orientation  that  were  acceptable  at  this  location.  The  orienta¬ 
tion  that  aj^ars  to  be  the  best  is: 

CX  -  90.  /3»  90. 

This  angle  matrix  is  listed  in  Table  VIII.  Note  that  modes 
11  and  12  have  completely  coupled.  Also,  while  angles  1-6  and  1- 

11  have  become  worse,  in  general,  the  coupling  characteristics 
have  improve.  Specifically,  observe  the  marked  improvement  be¬ 
tween  modes  6-11,  6-12,  7-11  and  7-12.  The  decoupling  has  also 
improved  throughout.  The  above  noted  improved  angles  are  circled; 
the  worse  angles  are  boxed. 

The  ORIENT  run  for  the  second  grouping  produced  no  orienta¬ 
tions.  More  exclusions  from  the  requirements  for  acceptance  need 
to  be  made  before  orientations  are  found.  The  new  sensor/actuator 
pair  picked  for  the  previous  grouping  is  not  acceptable  for  this 
grouping.  This  is  seen  from  the  fact  that  coupling  between  mode 

12  and  group  II  increases  instead  of  decreasing. 

The  next  objective  in  chemging  the  sensor /actuator  model  was 
to  create  new  groupings.  When  modes  11  and  12  were  decoupled, 
they  determined  two  of  the  three  groups.  When  the  modes  were 
coupled,  only  one  group  was  determined.  Two  of  the  possibilities 
for  groupings  that  were  found  (one  for  11  and  12  coupled,  one  for 
them  decoupled)  are  shown  in  Table  IX  and  Table  X.  The  mode 
groups  for  the  coupled  case  are: 

I:  1,  6,  7 
II:  2,  3,  4,  5 
III:  8,9,  10,  11,  12 


0 


The  mode  groups  £or  the  decoupled  case  are: 

I:  6,  1,  12 
II:  1,  5,  10,  11 
III:  2,  3,  4,  8,  9 

The  two  cases  studied  should  have  produced  most  of  the  pos¬ 
sible  mode  grouping  combinations;  however,  no  radical  improvements 
were  noted  with  any  sensor /actuator  model. 


If 


1 


0.000 


1 

0.000 

89.9870 

89.9990 

90.0098 

90.0541 

2 

89.9870 

0.0000 

67.4653 

175.3876 

116.8520 

3 

89.9990 

67.4653 

0.0000 

115.3141 

104.0988 

4 

90.0098 

175.3876 

115.3141 

0.0000 

60.6688 

5 

90.0541 

116.8520 

104.0098 

60.6688 

0.0000 

6 

106. 546  4| 

90.0434 

89.9875 

89.9621 

89.6876 

7 

79.1801 

89.9845 

90.0128 

90.0123 

89.9137 

8 

90.0119 

106.5687 

116.9916 

73.1913 

102.7659 

9 

90.0215 

116.2323 

91.5675 

66.4898 

68.5696 

10 

89.9951 

90.1561 

80.1261 

90.6814 

111.8731 

11 

107.6054| 

89.9918 

89.9969 

90.0092 

90.1359 

12 

107.7985 

89.9893 

89.9951 

90.0118 

90.1379 

1 

S. 

a 

Ifl 

11 

1 

79.1801 

90.0119 

90.0215 

89.9951 

|l07.6054| 

2 

89.9845 

106.5687 

116.2323 

90.1561 

89.9918 

3 

90.0128 

116.9916 

91.3675 

80.1281 

89.9969 

4 

90.0123 

73.1913 

66.4898 

90.6814 

90.0092 

5 

89.9137 

102.7659 

68.5696 

111.8731 

90.1359 

6 

50.7986 

90.0044 

89.9067 

90.1120  1 

^1.517^ 

7 

0.0000 

89.8663 

90.0804 

89.9508  1 

8 

89.8663 

0.0000 

93.8831 

38.0463 

90.1295 

9 

90.0804 

93.8831 

0.0000 

103.2217 

89.9727 

10 

89.9508 
- ^ 

38.0463 

103.2217 

0.0000 

90.0604 

11 

tf59.525ffl 

90.1295 

89.9727 

90.0604 

0.0000 

106.5464| 

90.0434 

89.9875 

89.9621 

89.6876 

0.0000 

50.7986 

90.0044 

89.9067 

90.1120 


12. 


12  0.59.830 


80.1314  89.9751  90.0614 


0.5245 


0.0000 


TABLE  IX 


Angle  Matrix  for  New  Grouping  fll  &  12  Coupled) 


12  1 


1 

0.0000 

91.3907 

91.7781 

2 

91.3907 

0.0000 

65.1063 

3 

91.7789 

65.1063 

0.0000 

4 

88.1414 

173.7748 

118.3582 

5 

91.4930 

112.3436 

99.4407 

6 

123.5268 

88.8425 

88.4525 

7 

109.2235 

87.8575 

87.3593 

8 

95.9688 

83.5529 

83.4629 

9 

93.3296 

104.5331 

82.8228 

.0 

89.4517 

91.2650 

81.9931 

.1 

83.3505 

102.7891 

106.3187 

.2 

83.5714 

102.7831 

106.3152 

1 

1 

i 

1 

108.2235 

95.9688 

93.3296 

2 

87.8575 

83.5529 

104.5331 

3 

87.3593 

83.4629 

82.8228 

4 

92.8263 

100.4695 

80.2084 

5 

87.8982 

82.2686 

65.1712 

6 

60.7614 

84.8955 

87.0308 

7 

0.0000 

80.9959 

85.1309 

8 

80.9959 

0.0000 

61.8256 

9 

85.1309 

61.8256 

0.0000 

.0 

90.7925 

79.1714 

103.8861 

.1 

99.4980 

159.8746 

121.6093 

.2 

99.9088 

159.8771 

121.6123 

A 

5. 

6. 

88.1414 

81.4930 

123.5268 

173.7748 

112.3636 

88.8425 

118.3582 

99.4407 

88.4525 

0.0000 

66.9967 

91.5642 

66.9967 

0.0000 

88.4894 

91.5642 

88.4894 

0.0000 

92.8293 

87.8982 

60.7614 

100.4695 

82.2686 

84.8955 

80.2084 

65.1712 

87.0308 

89.1670 

112.4157 

90.6307 

72.8646 

102.8915 

95.5664 

72.8712 

102.9008 

95.3021 

ID. 

11 

11 

89.4517 

.3505 

83.5714 

91.2650 

102.7891 

102.7831 

81.9931 

106.3187 

106.3152 

89.1678 

72.8646 

72.9712 

112.4157 

102.8915 

102.9009 

90.6307 

95.5664 

95.3021 

90.7925 

99.4980 

99.3388 

79.1714 

159.8746 

159.8771 

103.8861 

121.6093 

121.6123 

0.0000 

84.9988 

85.0043 

84.9988 

0.0000 

^^0^68^ 

85.0047 

Co^68^ 

0.0000 

TABLE.  X 


Anale 

Matrix  for 

New  Grouoina  (11  & 

12  Decouoled) 

1 

2 

1 

1 

2 

S 

1 

0.0000 

90.8174 

90.7932 

88.9493 

122.0821 

2 

90.8174 

0.0000 

66.7463 

174.9228 

116.6676 

86.3539 

3 

90.7932 

66.7483 

0.0000 

116.1557 

105.8840 

86.4521 

4 

88.9493 

174.9228 

116.1557 

0.0000 

60.1485 

94.6762 

5 

86.7474 

116.6676 

105.8840 

60.0485 

0.0000 

104.7235 

6 

122.0821 

86.3539 

86.4521 

94.6762 

104.7235 

0.0000 

7 

108.3965 

87.8155 

87.9756 

92.7696 

98.6674 

57.6279 

8 

90.3337 

105.9797 

116.2828 

74.0025 

102.1147 

88.7901 

9 

89.8009 

116.2560 

91.8262 

66.5076 

71.1859 

90.8291 

10 

90.4521 

89.5529 

79.6591 

91.4364 

110.4766 

88.1533 

11 

39.4739 

87.9660 

89.4698 

92.3365 

95.2409 

103.3852 

12 

82.9303 

86.0350 

87.3952 

94.7842 

102.5730 

68.8043 

2 

E 

2 

IS. 

11 

12. 

1 

108.3965 

90.3337 

89.8009 

90.4521 

39.4739 

82.8303 

2 

87.8155 

105.9797 

116.2560 

89.5529 

87.9660 

86.0350 

3 

87.9759 

116.2828 

91.8262 

79.6591 

89.4698 

87.3952 

4 

92.7696 

74.0028 

66.5076 

91.4364 

92.3365 

94.7842 

5 

98.6674 

102.1147 

71.1859 

110.4766 

95.2409 

102.5780 

6 

57.6279 

88.7901 

90.8291 

88.1533 

103.3852 

68.8043 

7 

0.0000 

89.1711 

90.6934 

88.8029 

81.0230 

120.3491 

8 

89.1711 

0.0000 

93.9749 

38.0012 

93.5142 

92.3150 

9 

90.6934 

93.9749 

0.0000 

103.3329 

90.4969 

91.2274 

10 

88.8029 

38.0012 

103.3329 

0.0000 

94.1676 

92.3073 

11 

81.0230 

83.5142 

90.4969 

94.1676 

0.0000 

(89.9243 

12 

120.3791 

92.3150 

91.2274 

92.3073 

(8^9243) 

0.0000 

The  results  produced  by  changing  the  original  model  orienta¬ 
tions  are  listed  in  Table  XI.  The  table  lists  the  new  angles 
produced  and  the  delta  change  from  the  old  angles.  The  plus  and 
minus  signs  indicate  favorable  or  unfavorable  chemges  respectively 
(i.e.r  minus  sign  if  new  angle  is  now  more  coupled,  but  the  desire 
was  to  be  more  decoupled).  The  results  show  six  of  nine  favorable 
changes  overall;  however,  there  are  five  of  six  favorable  changes 
when  the  desire  was  to  couple  (first  three  and  last  three).  It 
seems  to  be  easier  (I  don't  know  why)  to  couple  the  angles  than  to 
decouple  the  angles. 

The  magnitude  of  the  changes  are  small  especially  for  the 
three  cases  where  decoupling  was  the  object.  Yet,  the  largest 
changes  (8.5°  and  8.1°)  are  significauit  angle  movements  when  it  is 
remembered  that  only  one  sensor/actuator  orientation  was  ch2mged 
for  each  case.  Note  that  when  the  angle  changes  are  unfavorable, 
the  magnitude  of  the  changes  are  small  (0.3  is  the  largest). 


Mms  Model 


The  angle  matrices  produced  from  the  mass  model  changes  are 
cataloged  in  Table  XII  through  Table  XVI.  Note  that  while  some 
changes  in  angles  did  occur,  the  original  groupings  stayed  the 
same.  This  can  be  seen  by  comparing  Tables  IV  with  these  new 
t2d>lee. 

The  greatest  angle  changes  were  found  in  Table  XVI  (angles 
are  circled).  A  10^  change  occurred  between  modes  9-10.  Modes  2- 
10  and  4-10  had  6^  changes.  All  other  angles  in  the  model  changes 
had  angle  changes  of  less  than  3^.  These  results  indicate  that 
mode  10  is  more  sensitive  than  the  other  modes  to  the  tertiary 
mirror  mass  model  changes  that  caused  Table  XVI.  That  model 
change  was  the  one  with  equal  masses  at  the  corners  of  the  ter¬ 
tiary  mirror.  However,  in  general,  the  structural  model  is  insen¬ 
sitive  with  respect  to  modal  angles  to  the  tertiary  mirror  lumped 
mass  model  changes  tried. 
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i 

i  Q 
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IftBLB  m 

Matcix  with  2%  Hass. 

Change 

1 

2 

1 

A 

5 

6. 

1 

0.0000 

89.9877 

89.9990 

90.0086 

90.0890 

56.7100 

2 

89.9877 

0.0000 

67.2166 

175.3286 

116.8176 

89.9474 

3 

89.9990 

67.2166 

0.0000 

115.6233 

104.0412 

90.0152 

4 

90.0086 

175.3289 

115.6233 

0.0000 

60.6811 

90.0465 

5 

90.0880 

116.8176 

104.0412 

60.6811 

0.0000 

90.3062 

6 

55.7100 

89.9474 

90.0142 

90.0465 

90.3062 

0.0000 

7 

72.4413 

90.0340 

89.9624 

89.9781 

89.9860 

61.0909 

8 

90.0469 

106.7655 

116.8679 

73.0010 

102.7656 

89.9273 

9 

90.0131 

116.3286 

91.6029 

66.4258 

68.6255 

90.1274 

10 

90.0053 

90.0124 

80.1219 

90.8335 

112.0084 

89.8454 

11 

39.3911 

87.8653 

89.4109 

92.4748 

96.6414 

75.1253 

12 

81.3855 

87.9198 

89.2806 

92.3748 

95.6397 

107.0160 

2 

£ 

1 

11 

12 

"1 

1 

72.4413 

90.0469 

90.0131 

90.0053 

39.3911 

81.3855 

•  1 

2 

90.0340 

106.7655 

116.3296 

90.0124 

87.8653 

87.9198 

-  .. 

3 

89.9624 

116.8679 

91.6029 

80.1219 

89.4109 

89.2806 

4 

89.9781 

73.0010 

66.4258 

90.8335 

92.4749 

92.3748 

5 

89.9860 

102.7656 

68.6255 

112.0084 

96.6414 

95.6397 

— 

6 

61.0909 

89.9273 

90.1274 

89.8454 

75.1253 

107.0160 

7 

0.0000 

90.0882 

89.8507 

90.0423 

99.278 

52.6260 

8 

90.0882 

0.0000 

93.6552 

38.0181 

93.5231 

93.1494 

— 

9 

89.8507 

93.6552 

0.0000 

103.4248 

90.5090 

90.7426 

10 

90.0423 

38.0181 

103.4248 

0.0000 

94.1647 

93.5417 

11 

99.5278 

93.5231 

90.5090 

94.1647 

0.0000 

89.9618 

-  - 

12 

52.6260 

93.1494 

90.7426 

93.5417 

89.9618 

0.0000 
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Anale 

IftBLE  XIII 

Matrix  with  10%  Mass 

Change 

1 

2 

1 

i 

2 

6. 

1 

0.0000 

89.9880 

89.9989 

90.0081 

90.0912 

123.6228 

2 

89.9880 

0.0000 

66.2873 

175.1010 

116.7026 

90.0511 

3 

89.9989 

66.2873 

0.0000 

116.7874 

104.1296 

89.9839 

4 

90.0081 

175.1010 

116.7874 

0.0000 

60.7141 

89.9552 

5 

90.0912 

116.7026 

104.1296 

60.7141 

0.0000 

89.6929 

6 

123.6228 

90.0511 

89.9839 

89.9552 

89.6929 

0.0000 

7 

72.1172 

90.0348 

89.9626 

89.9777 

89.9922 

119.4366 

8 

90.0438 

107.4692 

116.4368 

72.3223 

102.7904 

90.0741 

9 

90.0143 

116.6894 

91.7448 

66.1947 

68.8212 

89.8733 

10 

90.0006 

89.5301 

80.1183 

91.3542 

112.4773 

90.1582 

11 

40.2559 

87.7743 

89.4148 

92.5773 

96.7645 

104.7872 

12 

82.4133 

87.8332 

89.2694 

92.4708 

95.7364 

71.8194 

2 

S. 

2 

12 

U 

12 

1 

72.1172 

90.0438 

90.0143 

90.0006 

40.2559 

82.4133 

2 

90.0348 

107.4692 

116.6894 

89.5301 

87.7743 

87.8332 

3 

89.9626 

116.4368 

91.7448 

80.1183 

89.4148 

89.2694 

4 

89.9777 

72.3223 

66 .1947 

91.3542 

92.5773 

92.4708 

5 

89.9922 

102.7904 

68.8212 

112.4773 

96.7645 

95.7364 

6 

119.4366 

90.0741 

89.8733 

90.1582 

104.7872 

71.8194 

7 

0.0000 

90.0825 

89.8513 

90.0394 

99.8400 

54.3089 

8 

90.0825 

0.0000 

92.8640 

37.9341 

93.6165 

93.2392 

9 

89.8513 

92.8640 

0.0000 

104.0881 

90.5053 

90.7533 

10 

90.0394 

37.9341 

104.0881 

0.0000 

94.2277 

93.5884 

11 

99.8400 

93.6165 

90.5053 

94.2277 

0.0000 

90.4336 

12 

54.3089 

93.2392 

90.7533 

93.5884 

90.4936 

0.0000 
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TABLE  XIV 


a 


Aggie  MatilA.  xith  Mass  Change 


1 

2 

1 

A 

5. 

6. 

1 

0.0000 

89.9883 

89.9986 

90.0076 

90.8951 

124.0894 

2 

89.9863 

0.0000 

65.2012 

174.8062 

116.6249 

90.0492 

3 

89.9988 

65.2012 

0.0000 

118.1720 

104.2001 

89.5824 

4 

90.0076 

174.8065 

118.1720 

0.0000 

60.6928 

89.9576 

5 

90.0951 

116.6249 

104.2001 

60.6928 

0.0000 

89.6903 

6 

124.0894 

90.0492 

89.9824 

89.9576 

89.6903 

0.0000 

7 

71.6635 

90.0359 

89.9628 

89.9772 

90.0010 

120.1351 

8 

90.0400 

108.2412 

115.9570 

71.5885 

102.8423 

90.0760 

9 

90.0156 

117.1463 

91.9086 

65.9119 

69.0926 

89.8745 

10 

89.9948 

88.8619 

80.1204 

92.8619 

113.0416 

90.1624 

11 

139.6665 

92.3373 

90.5802 

87.2966 

83.0830 

75.3050 

12 

83.6862 

87.7259 

89.2556 

92.5902 

95.8572 

70.4391 

2 

a 

& 

in 

11 

12 

1 

71.6635 

90.0400 

90.0156 

89.9949 

139.6665 

83.6862 

2 

90.0365 

108.2412 

117.1463 

88.8619 

92.3373 

87.7259 

3 

89.9628 

115.9570 

91.9096 

80.1204 

90.5802 

89.2556 

4 

89.9772 

71.5885 

65.9119 

92.0619 

87.2966 

92.5902 

5 

90.0010 

102.8423 

69.0926 

113.0416 

83.0830 

95.8572 

6 

120.1351 

90.0760 

89.8745 

90.1624 

75.3050 

70.4391 

7 

0.0000 

90.0755 

89.8518 

90.0362 

79.8297 

56.4074 

8 

90.0755 

0.0000 

91.8593 

37.9764 

86.2731 

93.3463 

9 

89.8518 

91.8593 

0.0000 

104.9249 

89.5103 

90.7575 

10 

90.0362 

37.9764 

104.9249 

0.0000 

85.6344 

93.6457 

11 

79.8297 

86.2731 

89.5103 

85.6944 

0.0000 

88.9311 

12 

56.4074 

93.3463 

90.7575 

93.6457 

88.9311 

0.0000 
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TABLE  XV 


1 

2 

1 

1 

5. 

1 

0.0000 

89.9876 

89.9991 

90.0088 

90.0885 

56.6687 

2 

89.9876 

0.0000 

67.9385 

175.4363 

116.9616 

89.9447 

3 

89.9991 

67.9385 

0.0000 

114.7622 

103.8782 

90.0138 

4 

90.0088 

175.4363 

114.7622 

0.0000 

60.5714 

90.0495 

5 

90.0885 

116.9616 

103.8782 

60.5714 

0.0000 

90.3092 

6 

56.6687 

99.9447 

90.1328 

90.0495 

90.3092 

0.0000 

7 

72.3954 

90.0351 

89.9620 

89.9773 

89.9845 

61.1721 

8 

90.1492 

106.3473 

116.6995 

73.4728 

102.7894 

89.9300 

9 

90.0129 

116.3590 

91.3043 

66.3991 

69.1853 

90.1267 

10 

90.0061 

88.2308 

79.9539 

92.3968 

112.6655 

89.8403 

11 

140.9311 

91.9666 

90.5957 

87.6792 

83.4320 

104.8729 

12 

80.6971 

88.0763 

89.2928 

92.2303 

95.5756 

105.9443 

2 

S. 

2 

IQ. 

11 

12 

1 

72.3954 

90.0492 

90.0129 

90.0061 

140.9311 

80.6971 

2 

90.0351 

106.3473 

116.3590 

88.2308 

91.9666 

88.0763 

3 

89.9620 

116.6995 

91.3043 

79.9539 

90.5957 

89.2928 

4 

89.9773 

73.4728 

66 .3990 

92.3958 

87.6792 

92.2303 

5 

89.9845 

102.7894 

69.1853 

112.6655 

83.4320 

95.5756 

6 

61.1721 

89.9300 

90.1267 

89.8403 

104.8729 

105.9443 

7 

O.OOOU 

90.0886 

89.8484 

90.0523 

80.7988 

51.4492 

8 

90.0886 

0.0000 

92.3421 

39.1026 

86.5713 

93.0572 

9 

89.8484 

92.3421 

0.0000 

105.7099 

89.7266 

90.5322 

10 

90.0523 

39.1026 

105.7099 

0.0000 

85.9199 

93.4492 

11 

80.7988 

86.5713 

89.7266 

85.9199 

0.0000 

90.4827 

12  51.4492  93.0572  90.5322  93.4492  90.4827 


0.0000 


TABLE  XVI 


i 


I 

y'- 
\  • 

i 

i 

pm 

i 

i 

I 

i 


a 


Angle  Matrix  with  Equal  Mass  Corners 


1 

2 

1 

A 

5. 

£ 

1 

0.0000 

89.9880 

89.9988 

90.0081 

90.0946 

123.2426 

2 

89.9880 

0.0000 

65.7677 

175.0560 

116.6448 

90.0466 

3 

69.9988 

65.7677 

0.0000 

117.3911 

104.4288 

89.9820 

4 

90.0081 

175.0560 

117.3911 

0.0000 

60.7695 

89.9596 

5 

90.0946 

116.6448 

104.4288 

60.7695 

0.0000 

89.6890 

6 

123.2426 

90.0466 

89.9820 

89.9596 

89.6890 

0.0000 

7 

72.3479 

90.0355 

89.9619 

89.9762 

89.9951 

119.5189 

8 

90.0471 

108.2480 

114.3825 

71.7586 

102.6075 

90.0740 

9 

90.0147 

116.6529 

91.6308 

66.2921 

69.4423 

89.8772 

10 

89.9984 

(84^34^ 

79.9181 

(96^01^ 

114.6969 

90.1752 

11 

40.4695 

88.0756 

89.3477 

92.3070 

96.6715 

104.6088 

12 

83.0063 

88.0505 

88.2420 

92.2379 

95.6539 

70.8076 

1 

2 

2 

12 

11 

12 

1 

72.3479 

90.0471 

90.8147 

89.9984 

40.4695 

83.0063 

2 

90.0355 

108.2480 

116.6529 

^4^34^ 

88.0756 

88.8905 

3 

89.9619 

114.3825 

91.6308 

79.6181 

89.3477 

89.2420 

4 

89.9762 

71.7586 

66.2921 

92.3070 

92.2379 

5 

89.9951 

102.6075 

69.4423 

114.6969 

96.6715 

95.6539 

6 

119.5189 

90.0740 

89.8772 

90.1752 

104.6089 

70.8976 

7 

0.0000 

90.0737 

89.8516 

90.0595 

100.1910 

55.3692 

8 

90.0737 

0.0000 

89.5405 

41.7858 

93.5619 

93.1830 

9 

89.8516 

89.5405 

0.0000 

^3.100 

90.0662 

90.3736 

10 

90.0595 

41.7858 

^13a^ 

0.0000 

94.0219 

83.3597 

11 

100.1910 

93.5619 

90.0662 

84.0219 

0.0000 

90.8776 

12 

55.3692 

93.1829 

90.3736 

93.3597 

90.8773 

0.0000 
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VI  Conclusions  and  ReccMmoendatlons 

By  separating  the  system  matrix  modes  into  orthogonal  groups, 
the  system  will  be  inherently  decoupled.  It  is  apparent  that  the 
judicious  choice  of  sensor/actuator  location  and  orientation  can 
achieve  acceptable  orthogonal  groups.  It  is  also  obvious  that  a 
single  misplaced  sensor/actuator  pair  can  destroy  a  previous  modal 
grouping  set.  The  practical  result  is  that  modal  groups  can  only 
be  partially  decoupled  because  of  the  difficulty  in  finding  the 
ideal  sensor /actuator  model. 

The  attempts  to  decouple  modes  11  and  12  show  the  ease  with 
which  a  single  coupling  angle  can  be  changed.  However,  the  diffi¬ 
culties  in  finding  an  orientation  that  will  change  only  one  euigle 
and  leave  the  others  unchanged  was  also  demonstrated,  ihe  problem 
of  trying  to  improve  all  coupling  characteristics  can  only  be 
solved  by  concentrating  on  improving  one  angle  at  a  time.  By  this 
method,  a  systematic  study  should  produce  the  desired  sensor/ 
actuator  model. 

The  studies  of  the  mass  model  changes  showed  little  effect  on 
the  modal  angles. 

The  next  step  in  the  study  of  sensor/actuator  placement 
should  be  that  of  more  modes.  In  this  study,  twelve  modes  are 
separated  into  groups.  The  methods  of  this  study  can  be  applied 
to  a  larger  number  of  modes.  However,  a  more  rigorous  method  of 
determining  the  modal  groupings  needs  to  be  made. 

A  possible  method  is  a  variation  of  the  method  that  ORIENT 
uses  (requirements  of  acceptance).  The  angle  matrix  could  be  in- 
puted  into  a  program  that  would  output  possible  groupings  based  on 
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some  given  criteria.  This  criteria  would  be  similar  to  that  used 
in  ORIENT;  all  angles  coupled  within  some  limit  would  be  in  the 
same  group  and  those  angles  decoupled  would  be  in  different 
groups.  The  angles  that  don't  fit  the  criteria  would  be  listed  as 
problem  angles  for  each  grouping.  These  angles  would  then  be 
subject  to  correction  using  ORIENT. 

I  feel  that  the  conditions  of  changing  one  angle  without 
changing  any  other  angle  needs  to  be  studied.  If  it  can  be  deter¬ 
mined  when  this  condition  is  applicable,  it  will  be  a  powerful 
tool.  This  method  could  be  applied  to  every  angle  until  all 
groups  were  perfect.  Implementation  of  this  method  would  be  in 
ORIENT'S  requirements  of  acceptance. 

One  other  step  needs  to  be  taken:  experimental  evaluation. 
The  case  in  controlling  a  structure,  when  an  added  sensor/actuator 
improves  the  modal  groups'  coupling,  can  be  studied  experimen¬ 
tally.  Thus,  once  a  sensor /actuator  model  is  found  to  be  accept¬ 
able,  it  can  be  tested  to  see  if  it  is  experimentally  acceptable. 
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The  listing  below  is  the  NASTRAM  Data  Deck  for  the  CSDL  II , 
Revision  3  large  space  structure.  Cards  of  interest  are:  BIGR, 
GRID,  CBAR,  C0NM2  and  PBAR. 

The  input  data  card  EIGR  stands  for  Real  Eigenvalue  Extrac¬ 
tion  Data.  The  method  used  is  the  inverse  power  (INV)  method  with 
symmetric  matrices.  The  frequency  range  of  the  eigenvalues  is 
given  with  the  next  two  numbers,  0-14hz.  The  estimate  of  the 
number  of  eigenvalues  in  the  frequency  range  is  given  next,  while 
the  following  number  asks  that  all  the  eigenvalues  be  given. 

GRID  is  the  input  card  which  defines  the  location  of  all  the 
nodes  of  the  model.  This  card  is  very  useful  in  determining  dif¬ 
ferences  between  similar  models. 

The  CBAR  card  is  Simple  Beam  Element  Connection.  This  card 
identifies  each  bar  element,  and  gives  the  nodes  that  are  con¬ 
nected  to  its  ends. 

OOSM2  is  the  card  that  was  varied  in  this  investigation.  It 
stands  for  Concentrated  Mass  Element  Connection.  This  card  allows 
lumped  masses  to  be  put  at  the  nodes.  The  mass  to  be  changed  is 
the  third  number  on  the  card.  On  some  of  the  cards,  the  moments 
of  Inertia  are  also  given  following  the  lumped  mass. 

Simple  Beam  Properties  are  listed  in  PBAR.  The  properties 
listed  are:  area,  moments  of  inertia,  and  the  torsional  constant. 
This  card  could  have  been  used  to  alter  any  of  the  properties  of 
the  mirror  support  trusses. 
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This  program  was  used  to  pick  off  the  critical  modes  from  the 
eigenvectors  obtained  from  NASTSAM.  The  input  is  read  from  TAPE 
5,  and  the  output  (critical  eigenvectors)  is  listed  on  TAPE  6. 

The  input  file  contains  three  pieces  of  information.  First, 
the  number  of  critical  modes  (N).  Second,  the  list  of  eigenvector 
numbers  (ID),  and  lastly,  the  38  eigenvalues  and  eigenvectors. 
The  output  from  NASTRAN  needs  to  be  modified  so  that  only  the 
eigenvectors  and  eigenvalues  are  present. 

The  eigenvectors  are  placed  into  a  three  dimensional  array, 
HAT  (59,  6,  38).  The  desired  eigenvectors  are  then  written  into 
TAPE  6  using  the  ID  vector  to  identify  the  critical  ones. 

Following  is  a  description  of  the  variables  used  in  SELECT: 

ID  (I)  s  The  eigenvector  numbers  of  the  critical 

modes. 

M  «  Index,  set  equal  to  ID  for  each  mode. 

NAT(I,J,K)  «  Three  dimoisional  array  of  eigenvectors. 

nODE  (K)  ■  List  of  eigenvector  numbers. 

>  Number  of  critical  modes. 
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This  program  uses  the  sensor/actuator  model  and  the  critical 
modes  to  calculate  the  system  matrices  (B  and  C)/  and  the  coupling 
angles  between  the  system  modes. 

The  program  first  reads  the  sensor/actuator  model.  The  num¬ 
ber  of  critical  modes  (N)/  actuators  (NA)r  and  sensors  (NS)  are 
read  as  integers.  Next  the  orientation  angles  of  the  actuators 
ALPHA,  BETA,  and  GAMMA  are  read.  The  identities  of  the  nodes 
where  the  actuators  are  placed  are  listed  in  NODE.  Note  that  this 
array  allows  multiple  actuators  at  one  node.  The  sensors  are 
similarly  read. 

The  critical  modes  are  read  next  into  PHI.  To  find  the  value 
corresponding  to  a  particular  node's  degree  of  freedom,  the  ID 
vector  is  used.  This  vector  orders  the  nodes  of  the  structure: 


ID  (1-19) =1-19 
ID (20-34) =26-40 
ID (35-50) =42-57 
ID(51)=100 


ID  (52) =910 
ID (53) =1001 
ID (54) =1002 
ID (55) =1003 


ID (56) =1004 
ID (57) =1112 
ID(58)=2830 
ID (59) =3233 


The  calculation  of  the  system  matrices  uses  the  formulas 
given  in  the  text  for  elements  of  B  and  C.  The  only  complication 
involves  picking  out  the  proper  element  of  the  eigenvector.  This 
is  accomplished  using  NODE  and  ID.  NODE  identifies  the  node  in 
question,  while  ID  translates  this  identification  to  determine 
which  of  the  59  nodes  it  is.  This  determination  allows  the  proper 
PHI  value  to  be  picked  out. 


69 


After  the  system  matrices  are  calculated,  it  is  a  simple 
matter  to  determine  the  emgles  between  the  modes.  The  dot  product 
formula  given  in  the  text  is  used  to  make  the  calculations.  The 
norms  squared  of  the  rows  (B)  and  columns  (C)  are  given  by  BSQ  and 
CSQ  respectively.  Dot  products  are  CB  and  CC.  The  angles  are  put 
into  matrices  ANGB  and  ANGC. 

Following  is  a  description  of  the  variables  used  in  ANGLE: 

ALPHA  (I)  =  0<'  ,  the  angle  between  +x  axis  and  the 

actuator . 

ALPHS  (I)  a  0(,  the  angle  between  +x  axis  and  the 
sensor. 

ANGB  (I,J)  a  The  output  matrix  of  emgles  between  the 
rows  of  B. 

ANGC  (I,J)  -  out£X]t  matrix  of  angles  between  the 

columns  of  C. 

B  (I,J)  =  System  matrix  which  is  the  ^  matrix 

times  the  direction  cosine  matrix  of  the 
actuators. 

BETA  (I)  =  0f  the  angle  between  +y  axis  and  the 

actuator . 

BETS  (I)  -  the  angle  between  -t-y  axis  and  the 

sensor . 

BSQ  (I)  -  The  sum  of  the  squares  of  each  element 

in  a  row  of  B. 

C  (I,J)  B  System  matrix  which  is  the  direction  , 

cosine  matrix  of  the  sensors  times  the^ 
matrix 

CB  «  Ihe  dot  product  of  two  rows  of  B. 

CC  s  The  dot  product  of  two  columns  of  C. 

CON  =  Conversion  factor  ,  Tt/lQO 

CSQ  (I)  B  The  sum  of  the  squares  of  each  element 

in  a  column  of  C. 

GAMMA  (I)  >■  the  angle  between  -i-z  axis  and  the 

actuator . 

GAMMS  (I)  "  the  angle  bestween  tz  axis  and  the 

sensor . 


ID  (I) 


N 


NA 

NODE  (I) 
NODS  (I) 
NS 


PHI  (I  J) 
PI 


=  The  node  that  corresponds  to  the  rank 
number. 

-  The  number  of  critical  modes. 

=  Ihe  number  of  actuators. 

=  The  list  of  actuator  nodes. 

=  The  list  of  sensor  nodes. 

-  The  number  of  sensors 

=  Input  matrix  of  critical  modes. 

=  TT 
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c'  -  IG’  AM  A^. GLF  (  t  'lTPUT  *T  APE'^  APL^) 

:  PH:  C354  >  tAL--HAC^‘  )  ,Ht  rA(=:s  ),GAM«A(S  ^  » 

niMEX'ilON  ALPHG<-.  >tfV‘TS<  ^  5),GAMv'  (  )  ,r.'D"(5-  )»'  ’'D' 

^)rM^^:^IC^  A\GH<3.  »  ^ C  )  ,  A'.'GC  (  3 :  »  3  3  )  t  RE  0  (5"'  >  »C Q  ( ~  ) 

□  :  ME',  lie?.  B  (  3  C  ,  E  :•  )  «C  (  5~  •  3  .  ) 

“  L  AO  (  '1 »  *  )  I,  »  f.A  <  fi  I 
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V 
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■ 
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E'UMBE*- 
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m 

p 1  Nj  r  * .  * 

NUMBEP 
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E  EAnC  T*  * 
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'.vDE 

ALt'HA 

HEO  A 

00  1  1=1 

.  A 

p  p  I  \ • 

« 

o  p  1 1\  r  3  0 

3  .1  tL(  OE( 

r  )  f  ALPHA  (  I  ), 

BETAC I ) ,GAMVA(  I  ) 

' EAO<E»* 

)  (ALFHSd 

» 

' rA0(5,  * 

)  (  I) 

»I=1 fNO) 

"tlAOCb** 

XGAMMSd  )il  =  l»?i  ) 

7  EAO( 5f  * 

)<^^DSCI) 

POI\' *,• 

n 

PHI\' *f " 

p  s  r  ’-i '  * , " 

m 

OR!  ti,TA 

ILL  OF 

0E.’.S''P"  " 

p  u  f.  ■  *  ♦  " 

m 

P  :  'g '  *  f  - 

t* 

.<.OE 

ALPHA 

BF^A 

00  Z  I=l 

P  w  I  f;  *  ,  • 

m 

PPIM  3  0'-  tl  OG( 

D.ALPH  (i)f 

HLfSd  ) 

fGAMK- '■  (  I) 

->  EAL)(5  f  A 

CC  )< <:  HI  ( 

I*d),I=l »3GA)»d=l »\ 

) 

’ tAO( F, . 

)  (IO(  I  )  ♦ 

=  1 t5  '  ) 

GA«'<A  • 


GAM'^  A  ' 


^  1=3.  lAlt'  2Gb’  C 

c -'?•)='■  1/1  ■  e.c 

O'  3C  J  =  1  V  A 

ALPHA  tJ»  =  ALPHA (J) *C  f 
HPT A( J)=BrTA< J)«C  \ 
GAMMA(  J)  =  gamma  (J)  ‘C 


0.  IE  K=lfEE 

1C  :F(I0(K).EQ.'  ■  DE(  J)  »G  T.  2C 
''c  c:t'^'?‘:..UE 
K=P*i^-5 
n:  31  :=if'' 

0  (  I  ♦  J)  =  PHI  (K*  :  )  -C  •  S(  ALPHA  (d»  )  ♦P  H I  (  !<♦  I  1 1  )  *  Ci  '  (  BE  ’  A  (  d  )  ) 
1  ♦PHI  (K^;  ♦!  )*Cr-:  <GAMMA<d)  ) 

3r  CCNTIiiUE 


*0( 
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00  60  1=1, NS 
ALPHSC  )  =  ALPHi(I  )  *CON 
BETS( I >=BETS( : )*C  CN 
G  AMM:.  <I  )  =  GAMM:;  (  I)  *CON 
00  4C  K=l,59 

40  rF(ID<Kl.£Q,f.CO$(  I)  )GO  Tvl  5C 
50  CONTINUE 
K=6*K-5 
00  60  J=l,f 

C (I ,J)=PH! {K,J)*C  :S<ALPHS (I ) ) ♦PH! (K  +  J  , J)*C: S(BE 
1  ♦PH:(K*-r,J)*CO£('.AFM  (1)1 

“^-O  CONTINUE 
00  70  I=l,‘. 

a  :  Q  ( I » =  G . : 

00  7C  K=1,NA 

a  n  =  0(i  ,K)**2.  :>b:  q(Ii 

70  CONTINUE 
m  9C  1=1 »N 
OG  93  d  =  i,f; 

08=0.'" 

0  0  80  K  =  1,'A 
‘•0  CB=B(  I,K)*B<J,K)>CB 

A  NGB(  r  ,  J)=ACCS<CB/Sas  r<8SQCI»  )/,>Q  T(H  ;  Q(U>>  ) 
ANGB<I,d)=ANGe<I, J)/C  N 
93  CONTINUE 

00  IQO  1  =  1, f' 

C  >Q<  I  )  =  0,0 

o:  ICO  K=1,NS 

CSa(I)=C(K,I)**2.^*C3Q(I» 

100  CONTINUE 

0  3  120  1  =  1  ,' 

0  3  12C  Jrl,'' 

cc=o.: 

03  110  k=i,n; 

110  CC=C(K,I)*C(K,J)*CC 

ANGCC I,d»  =ACC£<CC/SQRTCC£Q<I ) >/0Q9  TCC ■ Q( J))) 
ANGCd,  *J)=ANGr<  I,vJ)/C  ’* 

120  CONTINUE 
PRIfc  2  0C 
P  R I N"  2  0  3 

PRINT  21C ,< 1,1=1, ' A) 

DO  500  1=1, ^ 

500  PRIN'  2  2C  ,I  ,(H(  I  ,  J)  ♦J  =  l  ,•  A  ) 

PRIN'  20C 
PRIN'  2C4 

PRINT  21C  ,(I  ,,r  =  i ,  .) 

0  3  63  3  I  =  l,f.O 

600  PFIN-'-  22U  ,I  ,(C<1  ,  J)  ,J=l,  .) 

PRIN'  20'" 

P^IN’’  201 

PRINT  2i:,(i,:=i,  ) 

DC  130  1=1, f 

13C  PRIN’  22C,I,(ANGO(r,J>,J=l,N) 


PRINT  20G 
PRINT  202 

PRINT  21C*(I.:=1.*  ) 

0"  14G  I=ltN 

140  PRINT  220  ,I  »(ANGC  (I  tJ)*  J=1  »»:) 

200  FORMATIlHl) 

210  f:?RMAT(1X,5X.I7,11I10,/) 

220  FORNATdX  ./tib  tl2F10 .4,/,tX*12FlC  .^) 

201  F0RMAT(lXi2CX«24HANGLFS  BLTyFfN  ROWS  OF  R) 

202  F0RHAT(lXf2CX*27HaNGLES  BFTUFFN  COLUV'.O  CF  C) 

203  FORMAT! IX, 2OX.XH0  MATRIX* 

2C4  FORMAT! 1X,20X*“HC  MATRIX* 

3  00  F  0RMAT!1X ,1 5 , i 1 1 ♦FI 1 .2 .2F 10 .2 > 

400  FGRMA'^!lF3E15-6* 


Appendix  D 

ORIENT 

ORIENT  is  a  modified  ANGLE  program  used  to  calculate  the 
angle  matrix  for  many  different  actuator  models.  Two  major  dif¬ 
ferences  between  ANGLE  and  ORIENT  exist.  ORIENT  has  an  added 
actuator  that  varies  its  location  and  orientation.  It  also  has 
requirements  of  acceptance  that  give  a  measure  of  improvement. 

ORIENT  first  calculates  the  B  matrix  without  the  added  actu¬ 
ator.  This  part  of  the  B  matrix  remains  constant.  Only  a  new 
column  needs  to  be  added  with  the  new  actuator.  This  new  column 
is  calculated  for  every  node  point  and  every  possible  orientation 
(step  at  5^  increments). 

After  the  new  column  is  C2aculated,  AN6B  is  determined  as  in 
ANGLE.  However,  after  the  angle  is  calculated  for  each  mode 
pair,  it  is  then  tested.  If  any  angle  does  not  meet  the  require¬ 
ments  of  acceptance,  the  program  immediately  steps  to  the  next 
possible  orientation  for  the  actuator.  If  every  angle  meets  the 
requirements,  ANGB  is  printed  out  for  the  particular  orientation 
and  location  of  the  actuator  at  that  time. 

This  version  of  ORIENT  establishes  the  requirements  of  ac¬ 
ceptance  by  observing  the  coupling  characteristics  between  the 
modes  of  different  groups.  The  requirements  differ  depending  on 
what  groups  the  modes  are  in.  If  the  modes  are  in  the  same  group, 
they  are  required  to  be  coupled.  If  they  are  in  different  groups, 
the  angle  has  to  indicate  decoupling.  This  facet  of  the  program 
is  implemented  by  reading  each  mode  number  into  a  group  (Nl,  N2, 
or  N3)  at  the  beginning  of  ORIENT. 


This  mode  group  identification  allows  great  flexibility. 
First,  a  mode  grouping  can  be  tried  by  just  assigning  the  modes  to 
different  groups.  Second,  the  threshold  for  coupling  and  decou¬ 
pling  can  be  reassigned  to  any  desired  level.  This  is  the  method 
with  which  problem  angles  are  deedt. 

Following  is  a  description  of  the  variables  used  in  ORIENT: 

A  «  Possible  0< orientation,  calculated  at 

5®  increments. 

ALPHA  (I)  -  0(r  the  angle  between  +x  axis  and  the 

actuator . 

ANGB  (I,J)  -  The  output  matrix  of  angles  between  the 

rows  of  B. 

ARG  s  The  cosine  of  the  angle  between  the 

rows  of  B,  argument  of  arc  cosine  equa¬ 
tion. 

B  (I,J)  s  System  matrix  which  is  the  (!>  matrix 

times  the  direction  cosine  matrix  of 
the  actuators. 

BE  s  Possible  0  orientation,  calculated  at 

5®  increments. 

BETA  (I)  =  the  angle  between  +y  axis  and  the 

actuator . 

BSQ  (I)  s  The  sum  of  the  squares  of  each  element 

in  a  row  of  B. 

CB  B  The  dot  product  of  two  rows  of  B. 

CX)N  =  Conversion  factor , 'ft /1 80. 

G  «  Possible  ^  orientation,  calculated  from 

A  and  BE. 

G2  =  180®  -  G 

GAMMA  (I)  ,  the  angle  between  +z  axis  and  the 

actuator . 

11  s  Index  of  modes  in  group  1. 

12  «  Index  of  modes  in  group  2. 

13  B  Index  of  modes  in  group  3. 


ID  (I) 

N 

N1  (I) 
N2  (I) 
N3  (I) 
NA 
NAA 
ND 


B  The  node  that  corresponds  to  the  rank 
number. 

s  The  number  of  critical  modes, 
s  First  group  of  modes. 

=  Second  group  of  modes, 
a  fniird  group  of  modes, 
s  The  number  of  actuators. 

=  NA  +  1. 

s  Index,  set  equal  to  ID  for  each  node. 

=  the  list  of  actuator  nodes, 
a  Input  matrix  of  critical  modes. 


NODE  (I) 
PHI  (I,J) 


Pf-  CG- AM  .  SIf’JT(I'';PL)T,0UTPUT,TAPE5«TAPE&) 

DIMf;  PHI  (35 <*1 3 C  )tALPHA(59>*BETA(59),GAMMA(5v» 

ni  ►f/i  N  .  OF.Cov  )  ,  rO(55)  .M  (5)  *^2(5)  ,N3(5) 
niMF.\.5:0'S  Af.GH{3:  t3:)  *850(5?  )  *8(30  #59  ) 

-:5A0(5»*  )F«*.  A 
-LADCj**  )  (f  1  (I  )♦  i  =1  *4) 
tA0(5, *)(;:(!),:  =  !, 5) 

FAD(5»* )» : =1 *3) 


-LAD (5 
tA0(5 
LAD  (5 
P  5  I  \ 

P9i:.  * 
rcT.,7* 
P= IF’ * 
LAD  (5 
•  LAD (5 
LAO  (5 
rAD(^ 
h  I * 
pp  Ir:  ■  * 
p  p  I * 
PR  ’  * 
PP 

DO  1  I 
PRI.\  * 
c  r;,\ 

■EAD(D 

EAO(5 


.«‘BER  OF  CRITICAL  MODES  = 


KU^-REP  OF  ACTUATOPS  = 


* ) (AlPHA( I  »  » :  =  1 tPA) 
*)(BFA(I )*I=i»GA> 

*  )(GA*'MA(I),:  =  l,i\.A> 

*  )  (F  OLd  >  t  1  =  1  *P  A) 


ORIE^TA’•IO^  OF  ACTUATCRL 


ALPHA 


GAMMA 


pRIl.  50'’  .1  *1.  ODE  (I  )*ALPHA(  :>  tBE’’ A(  I)  »GAM*^A(:  ) 

■EAD(0«<»  .C  )((°HI(  I»J)*  1=1*354  )*J  =  1»N) 

EAO(5**  KIO(  I )  *  :  =1  *59J 
■;=3.14i502£;5EO 
CPN=r I/IR  C.Q 
',AA  =  -  A-»l 
0:'  31'  0=1* ’.A 
AlFHA(  0)  =ALPHA(d)  *C 
HrTA(J)=HL'A(J)*C:N 
«,AfRA(  «jl  rGAK'<A(  J)  •C  N 
0  1  ^  K=i*5- 

!F(  *n(K)  .ED,''  '0E(0>  )GC  TP  2  0 

^':^T:^alL 

K-9,*'<-5 

DO  3:  I=i*r. 

H(I ,J)=t  HI(K*I )*CCS(ALFHA(J) )♦FHI (K*l *  I ) *CCS ( BET A ( 0 ) ) 
•»PH:  (K*"  »I  )  *CP'  (GA  M^Ai  j)) 

:j.  O-Pi  KKrl*}-' 

■.D=in(K‘< } 
l;  -^9  [A  =  1»<7 
I  A  A  =  L  *  ■  A  -  5 

a=:aa*c  ' 

c  :  -  i  r  H=i .  d 
■  HB  =  "-*  ’M-5 
ftrrlPR.C  N 

IFC  AA  +  I  DR.M  .-t:  )  GO  TO  31 
0  =  9'''  .C  *C. 

Gv.  T'.  32 
C.  LT  '  MjL 
C:CA=CG- (A) 

F(C  EA.l, .)  CCS  A  =  -C'‘ '  A 
.F(  ■  \(Hb  ).L  dC'OAl  GC  'P 


t3  =  AC"S( ;  IN(HE.)*COS(A''  INCCOSA/"  Ir^)(BE)  )  )) 
G2=1-G«C  N-G 
:^2  CC\t!!»!UE 

0.  IG=1*2 

IF(rC.EQ.i)  GG  ■  0  3? 

:F(I  AA  +  i  BB.EQ.VC  )  GC  TC  =?9R 
G  ~G2 

35  CCNTI\UE 
K=KK 

OC  S  1=1*'. 

B(I  ♦■,AA)=PHr(K,!  >  *CCG(A>-*-PHI  (K*i  ♦!  )*C03  (B) 

1  ♦PH!  (K+Ttl  J’oCCS  (G) 

6  CUNTr\UE 
DC  7'  1  =  1,'': 

H':q<  :  )=o  .0 

DO  71  K=i,*jAA 

n3Q<I>=H{I,K>**2.0^B:Q<I) 

7C  CuGTI^L’E 

00  iC  1=1,0 
O'  O';  j=i,o 
CH=  .'  .C 

DO  0'  K=1,'AA 

i-ir  CH=B<I  ,K)*BCJ,K»+CB 

ASG=CH/CQRT(HSQ(: ) )/oQF  T(B^Q<J)) 

:F(A  G.LE.l.  .  EC:  )  GC  TO  PI 
A^  GOd  ,J  >=C. 

GO  ”  '  8  2 
CON'IMiE 

A  ,GB<I , J) =ACC3(C8/S0P  T<BSQ< I ) ) /C QR T( BS Q( J ) ) ) 

•  2  C  TNT  .  ^.iUE 

A'  G8d  ,J)  =A  'G8(I  ,0)/C^N 
O'.  101  11=1,4 


101 

’F<I .ED. N1 <: 1) ) 
O'  1C 2  :2=i,:' 

GO 

^  r 

V 

104 

102 

IFC.EQ.G'^CO)) 

GC 

'  0 

IC'^ 

0  103  "3=1,3 

iC5 

:F(i.E0.i<3(:3)) 

GC 

r» 

116 

1C  A 

n  114  1 1=1,4 

114 

IF( J.EQ. M (i: ) ) 

r,.'!'  ■  ..  Ill 

GO 

T  .-4 

lie 

1  » t, 

A  V 

0  1 15  :  ;?  =  1,  r 

ll" 

IF(  J.EQ.  '  2  (d  )  ) 

GO  '  11. 

GC 

1 

IK 

1  Of 

D  116  13=1,3 

1 1  f 

F(d.EQ.G3<:^)) 

GC 

'  1? 

lie 

G'  Ill 

:U  TFCA*  GR(  i  ,J»  .1“^.  2.0.0- .ANGB(I,J).GT. 1133.0)  GO  To 

G  0  3  '  G  ' 

ll  i  ^FCA  GB(  :,J).G'..  6.1  .AND. A'.GBCI,J)  .17.07, 0)  GO  TO  on 

GC  '  O'"  ■" 

3  cv.nt:nue 

r'f  !  *.  c!'  ' 

PPI'.  •  5G  .  ,A/C::\,[3E/C0N,G/CGIN,ND 

Pr'I'.'  20  1 

r-  r,-  j ' ,  -  2 1  c  ,  ( i  ,  I  =  i  ,  ) 

o:  i  ’!c  :  =1 

130  f-PIM  22- ,T,  <  A-,GB(I  ,J),J=1,\) 


.CC  F  (IHl  ) 

21C  FC->'AT(lX»5V,I7,lin  ,/) 

12  C  F:R?<Ar(lX,/».5,l£Fir  .4,/,5Xfl2FlC.4) 

:oi  f CKMA''(ixt3;'<»£«»HAKGLE:?  be:twfe\  POW?  CF  H) 

'CT  F  :5FAT(1X»  I‘M:iltF11.2.2H0.2) 

F.-kva’  (1P3E1 

"Gc  Fc.F  f'ATCl  Xt2u  MfcHALPHA=,F4*0 .6H  8ETA=tF4.C  #  7H  GAMMA=tF4.3, 
1  -H  AT  f, -of:  »!'=♦//) 

•  ->■ 

c:N:'':\ut 

D 
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